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Abstract 



O 

Oh! 

^ QCD corrections are evaluated to FCNC processes like B — > X s uu, K — > 7rz/z/, I? — > / + / _ or 
^1 ■ -KTl — > i.e. to processes mediated by effective operators containing neutrino currents 

>- . or axial leptonic currents. Such operators originate from W-hox and Z-penguin diagrams 
^ \ in the Standard Model. QCD corrections to them are given by two-loop diagrams. We 
\ confirm results for those diagrams which are already present in the literature. However, our 
analytical expressions for the Wilson coefficients disagree, due to a subtlety in regulating 
spurious IR divergences. The numerical effect of the disagreement is rather small. The 
size of the perturbative QCD corrections compared to the leading terms depends on the 
renormalization scheme used at the leading order. It varies from to around 15% for a 
reasonable class of schemes. The uncertainty originating from uncalculated higher-order 
(three-loop) QCD corrections is expected to be around 1% . 
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1. Introduction 

Flavour Changing Neutral Current (FCNC) processes are known to be an important 
source of information concerning the Standard Model (SM) parameters, as well as a window 
towards new physics. However, predictions for their amplitudes are often plagued with un- 
certainties due to low-energy QCD dynamics. Processes dominated by VT-box or Z-penguin 
diagrams which involve both quarks and leptons form an important class of exceptions from 
this rule. Low-energy dynamics can be then efficiently factorized, and the theoretical pre- 
dictions can be made very precise for a given set of fundamental SM parameters. 

Phenomenologically, the most interesting processes belonging to this class are B — > X s vv, 
Kl — ■> 7r°z/z/, K + — > 7i + uu and B — > [I]]— FCNC decays involving neutrinos in the 

final state are obviously receiving contributions only from VT-box or Z-penguin diagrams, 
because the heavy bosons mediate all the neutrino interactions. In the B — > case, 
photonic penguin contributions vanish due to the electromagnetic current conservation for 
the lepton pair, while double-photon intermediate states are numerically irrelevant. The 
analogous K — > decays are, on the contrary, dominated by the CKM-favoured double- 
photon contributions and, in consequence, much less theoretically clean. 




Figure 1: One- loop diagrams which give leading contributions to the considered processes in 
the Standard Model. The charged would-be Goldstone boson is denoted by G. 

The W-box and Z-penguin diagrams which give rise to all those processes are presented 
in fig. 1. Since external momenta in these diagrams are of order m# or nix, i.e. much smaller 
than My/, Mz or m t , one can describe the processes in question by introducing an effective 
theory which is obtained from the SM by decoupling the heavy electroweak bosons and the 
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top quark 



(leptons and light quarks) + 



AG F 
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Here, Q n stand for dimension > 4 effective interactions like 



Qi 




v t :v tb (s la p L b)(h a p L i), 



(2) 



Q 2 




v t :v tb (s la p L b)(^ a p L u), 



(3) 



where g 2 is the SU(2) gauge coupling constant, Vij are the CKM matrix elements, and 
Pl — (1 — 7s)/2. Numerical values of their Wilson coefficients 



are found by matching the full and the effective theory Green's functions. The matching is 
performed perturbatively in gauge couplings and in (external momenta) /M\y. 

The two operators Qi and Q2 written explicitly above are the ones which give dominant 
contributions to B — > and B —>■ X s uu, respectively Operators relevant for K + — > n + uu, 
Kl — > 7[°uu and K F — > are obtained from Qi and Q2 by simply replacing 6-quarks by 

ci-quarks in the operators themselves and in indices of the CKM factors. 

There are two main reasons for introducing the effective theory description. One of 
them is summing up large QCD logarithms like [a s ln(M^/m^)] fc from all orders of the 
SM perturbation series. This is achieved by applying renormalization group equations to the 
Wilson coefficients C n . The other reason is the fact, that different-looking Feynman diagrams 
can give rise to very similar effective operators whose nonperturbative matrix elements can 
be related to each other. Such relations are used e.g. for extracting nonperturbative matrix 
elements in K — > -kvv from the measurement of K + — > 7i°u e e + 0. 

In the processes under consideration here, renormalization group evolution of the Wilson 
coefficients plays only a minor role. The dominant contributions originate from operators 
with only a single (V-A) quark current, like those in eqns. @ and (|3|). Such a current 
is not renormalized in a mass-independent scheme, because it is conserved in the limit of 
vanishing quark masses. Consequently, anomalous dimensions of these operators vanish, 
and RGE running of their Wilson coefficients may arise only at higher orders in (light 
masses) 2 /M^. Such higher-order running is numerically relevant (though subdominant) 



Cn = Cf + ^ + ... 



(4) 
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for K + -> ir+vv and (K L -> /+/ ) 5hort Distance only.Q 

So long as the RGE running is negligible, and the low-energy matrix elements can be 
found reliably, the main calculational effort is required at the point of perturbative matching 
of the full SM and the effective theory. Since masses of the light particles can be set to zero 
when matching is performed, Wilson coefficients of Qi and Q2 are the same as the Wilson 
coefficients of their <i-quark analogs. If the Wilson coefficients of Qi and Q2 were found only 
at the leading order in QCD, i.e. from the purely electroweak one-loop diagrams in fig. 1, 
they would contain an inherent uncertainty of order 10% due to two-loop gluonic corrections. 
If those two-loop corrections were not calculated, this uncertainty would often become the 
main theoretical uncertainty in the final predictions for branching ratios, for given values of 
the SM parameters. 

A two-loop calculation of the next-to-leading matching conditions for Qi and Q2 was 
completed several years ago Jl|. It has not been checked by any other group so far. In the 
present paper, we recalculate all the next-to-leading contributions to the Wilson coefficients 
of Qi and Q2- We confirm the results of ref. JTJ] for all the necessary two-loop Feynman 
diagrams. However, our final expressions for the Wilson coefficients disagree, due to a 
subtlety in regulating spurious IR divergences. 

Our paper is organized as follows: In the next two sections, we present a detailed descrip- 
tion of evaluating Wilson coefficients of Q\ and Q2. Section 2 is devoted to the leading-order 
(one-loop) matching. The QCD corrections which arise at two loops are taken into account 
in section 3. Section 4 contains a discussion of numerical significance of the QCD corrections. 

2. Leading-order matching for Q\ and Q 2 . 

Leading-order contributions to the Wilson coefficients C\ and C2 can be found by requiring 
equality of perturbative off-shell amplitudes generated by the full Standard Model and the 
effective theory. We need to consider b — ► sl + l~ and b —> svv transitions in the cases of 
Qi and Q2, respectively. On the Standard Model side, we have to calculate the one-loop 
diagrams from fig. 1 for vanishing external momenta. On the effective theory side, one needs 
to include only tree-level diagrams with four-fermion vertices corresponding to either Q% or 
Q2. 

2 It arises from diagrams with double effective operator insertions and charm-quark loops. They corre- 
spond to W-boxes and Z-penguins with charm quarks on the full Standard Model side. Such diagrams are 
suppressed by m^/M^ but enhanced by CKM angles with respect to the leading top-quark contributions 
1- 
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For purposes of the next section, we shall perform the matching in D = 4 — 2e dimensions, 
i.e. without using any identity from only 4-dimensional Dirac algebra. In consequence, two 
other operators need to be introduced explicitly on the effective theory side 

Q2 = -^KVaisja^^PLb^u^^PLu) - 16Q 2 . (6) 

The superscript "E" stands for the name "evanescent" which originates from the fact that 
such operators vanish in 4 dimensions due to the identity 

7 75- (7) 

This identity cannot be analytically extended to D dimensions. Operators Qf and Qf arise 
in D-dimensional matching, because W-box diagrams in fig. 1 contain triple products of 
Dirac matrices which cannot be reduced to anything simpler in D-dimensions. 

Elementary calculation of the one-loop diagrams in fig. 1 with vanishing external mo- 
menta in the Feynman-'t Hooft gauge gives us the following results for the leading-order 
contributions to the considered Wilson coefficients^ 

Cf = C (x t )-B {x t )+O(e), (8) 

Cf = C {x t )-AB {x t ) + O{e), (9) 

Cf (0) = C^ = -\B (x t ) + O(e), (10) 

where x t = ntf/M^ and 

Bn(x) = — rrrlnx ; -, (11) 

ov ; 4(x - l) 2 4(.x - 1)' v ; 

. . 3x 2 + 2x . x 2 - 6x 

Co(x) = - —\nx + - -. 12 

6{X — l) z 8{x — 1) 

The functions Bq and Co originate from P^-box and Z-penguin diagrams, respectively. 



3. Next-to-leading order matching for Q\ and Q 2 . 

Before performing the NLO matching, we need to learn the structure of one-loop UV 
counterterms on the effective theory side. Evaluating the UV QCD counterterms can be done 

separately in the lepton (Qi,Qf) and neutrino {Q21Q2) sectors, thanks to their different 

3 Here, we ignore contributions suppressed by m^/Myy. They are relevant in K + — > tt + vv, i.e. for the 
d-quark analog of Q2, because of their CKM enhancement. 
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leptonic content. In each of the two sectors, bare quantities are replaced by the QCD- 
renormalized ones as follows: 

CQ + C E Q E -> (CZ NN Q + CZ NE Q E + C E Z EN Q + C E Z EE Q E ) . (13) 

Here, Z^ denotes the usual quark wave-function renormalization constant. The remaining 
renormalization constants are found in the MS scheme from the following two conditions 
which the n-loop effective theory amplitudes have to satisfy [ffOj: 



Renormalized amplitudes proportional to the coefficient C of the "normal" operator 
Q have to be finite in the limit D — ► 4. Counterterms which make them finite can 
contain nothing but l/e k poles, with 1 < k < n. 

Renormalized amplitudes proportional to the coefficient C E of the "evanescent" oper- 
ator Q E have to vanish in the limit D — ► 4. Counterterms which make them vanish 
can contain nothing but l/e k poles, with < k < n — lif they are proportional to the 
"normal" operator, and with 1 < k < n if they are proportional to the "evanescent" 
operator. This means that uniquely defined finite counterterms occur in this case, too. 




Figure 2: One-loop diagram which determines renormalization constants in the effective 
theory. 

Recovering the one-loop renormalization constants in the present case is extremely easy, 
because there is only a single one- loop diagram to consider. It is shown in fig. 2. The 
four-fermion vertex in this diagram is either the "normal" operator Q or the "evanescent" 
operator Q E . Since at this point we are interested only in the 1/e UV poles originating from 
loop momentum integrals, we set the external momenta and the physical masses to zero, but 



keep a common mass parameter as the IR regulator in all the propagator denominators (TT 
We find 



Z NN = 1, Z NE = 0, Z EN = -^r + 0(a 2 s ), Z EE = l + 0(a 2 s ), (14) 
with r = 32 and r = —32 in the lepton (Qi, Qf ) and neutrino (Q2, Q E ) sectors, respectively. 



The above equations mean that the renormalization constant was always enough to 
remove the 1/e pole from the one-loop diagram in fig. 2. Actually, the first two equations 
are true to all orders in QCD, because the quark current in Q is conserved for vanishing 
quark masses, as already mentioned in the introduction. The important part of eqn. (|n|) is 
the finite renormalization constant Zen which is going to affect our final results. 




Figure 3: Two- loop QCD corrections to the W-hox. diagrams from fig. 1. 




Figure 4: Two- loop QCD corrections to the Z-penguin diagrams from fig. 1. 




Figure 5: Diagrams with quark mass counterterms. 

We are now ready to perform the NLO matching, i.e. to find C{ and C% . On the full 
Standard Model side, one needs to calculate the two-loop diagrams shown in figs. 3 and 4, 
as well as the diagrams with quark mass counterterms in fig. 5. We ignore the quark wave- 
function renormalization now. This is allowed so long as one does it on the effective theory 
side of the matching equation, too. As in the previous section, it is enough to calculate all 
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the diagrams in the limit of vanishing external momenta. The calculation is particularly 
simple when masses of all the light particles are set to zero. Then, the 1PI parts of the 
diagrams depend on masses only via the ratio Xt, and they are very easy to handle with 
recurrence relations for two- loop vacuum integrals 0. 

Setting light masses to zero may lead to generation of spurious IR divergences. They 
should cancel with similar divergences on the effective theory side, so that the final results 
for the Wilson coefficients are finite in the limit D — > 4. This cancellation occurs only 
after comparing the full and effective theory amplitudes. Thus, one needs to calculate both 
amplitudes in D dimensions, even after their UV renormalization. This was the reason why 
we needed to perform the matching for Qf in the previous section. 

The sum of the diagrams in figs. 3, 4 and 5 turns out to be finite in the limit D — > 4. 
However, as we shall see in the following, it contains finite terms originating from multipli- 
cation of 1/e spurious IR divergences with 0(e) terms from the Dirac algebra. 



Calculating the effective theory side of the next-to-leading matching condition is partic- 
ularly simple, because we have set all the light masses to zero on the full theory side. Thus, 
we now need to do the same on the effective theory side, where only light particles occur. 
Consequently, all the loop diagrams on the effective theory side vanish in dimensional regu- 
larization. What remains are only tree-level diagrams with insertions of either the original 
operator vertices or the UV counterterms. The matching equation takes then the following 
form: 

= ( const -) >< { [ c(0) + £ ( c{1) + rcm )} (Q)»~ + <<A~} • 

(15) 

We have now consistently ignored the quark wave-function renormalization on both the full 
and effective theory sides. 

Vanishing of all the loop diagrams on the effective theory side means, in particular, that 
all the UV divergences in those diagrams cancel with spurious IR divergences (which arise 
due to setting all the external momenta and masses to zero). Thus, UV counterterms on 
the effective theory side actually reproduce the spurious IR divergences. The spurious IR 
divergences are the same on the full and effective theory sides, and cancel in the matching 
equation. In our case, there are no 1/e poles involved in this cancellation. We only have 
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the finite "rC E( - QS) " term. However, 1/e poles would occur in a generic case.0 Thus, the 
"rC^ )" term on the effective theory side can be interpreted as the one which cancels finite 
terms originating from multiplication of 1/e spurious IR divergences with 0(e) terms from 
the Dirac algebra on the full theory side. 

Since we do not need to recover the NLO Wilson coefficients of the evanescent operators, 
we can take the limit D — > 4 in the matching equation (0). Then {Q E )tree vanishes, and we 
find 

cf = Ci(x t ) - B x {x u -1/2) + 0{e), (16) 
Cf = d(x t ) - m(x t , +1/2) + 0(e), (17) 



where 



_ , . x 3 + Ax , x 4 — x 3 + 20x 2 2 —3a; 4 — 3a; 3 — 35a; 2 + x 

Cl{x) = (^iy^ 2(1 - x)+ 2(x-i)3 lnx + 3(^1? lnx 

Ax 3 + 7x 2 + 29x dC (x) , u 2 , , 

+ ^ tt^ + 8a; " v y In 



3(x - l) 2 ' ~~ 9a; ~* M, 2 , ' 



-1/2) = —rLi 2 {l-x) + - — rlnx + — rr^-lna; 

[x — iy (x — iy 3(x — xy 

~3x 2 + 19x io aBo(x) /i 2 
3(a; — l) 2 ox 



Si (x, +1/2) = — L2 2 (l-a;) + - — - In x + — r^-hia; 

[x — iy (x — iy Q(x — i 



3a; 2 + 29a; aB (x) , u 2 



Here, /i stands for the renormalization scale at which the top quark mass is MS-renormalized. 
The functions Bi(x t , ±1/2) and Ci(x t ) originate from gluonic corrections to H^-boxes and Z- 
penguins, respectively. The difference between B 1 (x, —1/2) and Bi(x, +1/2) is proportional 
to Bq(x) 

B x {x, -1/2) - B x {x, +1/2) = 6B (x). (21) 

Our result for the function C\(x) is in perfect agreement with ref. |l|. However, the results 
for B\ (x, ±1/2) disagree. The disagreement is due to the term "rC E ^" in our matching 
equation flT5|). 

Since our results differ from the previously published ones due to a subtlety in regulating 

spurious IR divergences, it is worthwhile to cross-check them using light quark masses as 

4 They do occur e.g. when two-loop off-shell matching for 6 — > S7 is performed for vanishing light particle 
masses. 
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IR regulators. Calculation of the 2-loop SM diagrams is then somewhat more complicated. 
One-loop diagrams on the effective theory side are nonvanishing in this case, and one needs 
to calculate their finite parts. However, the evanescent operators are unimportant. The 
matching equation can be written in 4 dimensions. Both the effective and the full theory 
amplitudes depend on the light quark masses in the same way. Terms dependent on the 
light quark masses cancel out in the matching equation. The Wilson coefficients we obtain 
in the end are exactly the same as the ones already found in eqns. (|i~6|)-(f20|). Thus, our 
NLO matching results are cross-checked by calculating them twice with use of two different 
regulators for spurious IR divergences: dimensional regularization and light quark masses. 



4. Numerical significance of the NLO corrections 

Branching ratios of the decays B — > X s vv and Ki — > 7i uu are proportional to the 
squared Wilson coefficient | C 2 1 2 , while B — > l + l~ is proportional to | C 1 1 2 . Thus, numer- 
ical significance of the NLO matching in these decays can be seen in figs. 6 and 7 where 
| Cfc 1 2 are plotted as functions of the matching scale //. Dependence on this scale enters via 



x t = [mf IS (n)/Mw} 2 } 0£ s (fi) and explicitly in C^ 1 . We use 



nr 



pole 



175 GeV and a s {M z ) 



0.118. Solid lines present the full NLO predictions, while dashed lines correspond to the LO 
results. Dotted lines indicate what the NLO results would be if the "rC^ -*" term was not 
included in the matching equation. 



I°i|lo 



2 

NLO 



[i [GeV] 

Figure 6: \C\(n)\ 2 at the leading and next-to-leading order. 

One can see that the NLO QCD corrections to \Ck\ 2 range from to around 15% for 
a reasonable class of renormalization schemes, i.e. for /i e [Myy,m t ]. The remaining \i- 
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H [GeV] 

Figure 7: |C*2(/x) | 2 at the leading and next-to- leading order. 

dependence at NLO is an estimate of uncalculated higher order (three-loop) QCD corrections. 
It is around 1%. The numerical effect of the "rC^ -*" term is seen to be rather small (around 
2%). To some extent, it is due to the fact that it originates from H^-boxes. Contributions 
from W-box.es tend to a constant at large x t , while Z-penguin contributions grow with x t , 
and are dominant for the physical value of x t ~ 5. Thus, from the numerical point of 
view, our results are a confirmation of what has been found in ref. QXJ and in the following 
phenomenological analyses (12 . 



The decays K + — > n + uu and Kl — > are somewhat more complicated because their 

amplitudes receive other significant contributions which are not proportional to the Wilson 
coefficients discussed in sections 2 and 3. Both of them are affected by 0{m 2 c / M^) terms 
which are CKM-enhanced with respect to the leading top-quark contributions. Those charm- 
quark loops are numerically as large as 20%-40% of the top-quark ones 0. In addition, the 
Kl — ■> mode is actually dominated by contributions from double-photon intermediate 

states which arise at higher order in a em , but involve u-quark loops and, in consequence, 
are strongly CKM enhanced with respect to the W-box and Z-penguin diagrams. Thus, 
QCD corrections to the top-quark loops are of mild importance for K + — > 7t + ujj, and totally 
unimportant in the Kl — > case. 



5. Summary 

We have evaluated QCD corrections to processes mediated by Z-penguin and W^-box 
diagrams. Two-loop matching of the full Standard Model and the effective theory was per- 
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formed off-shell, for vanishing external momenta and light particle masses. In consequence, 
calculation of the necessary two-loop Feynman diagrams was very easy, but the matching had 
to be performed in D dimensions, due to dimensionally regulated spurious IR divergences. In 
effect, Wilson coefficients of evanescent operators had to be found, and they affected the final 
physical results. Such a phenomenon is expected to take place in a generic process of heavy 
particle decoupling, when light particle masses are set to zero in the matching procedure. 

We have confirmed the results of ref. fl[] for the two-loop diagrams, but our final ex- 
pressions for the NLO Wilson coefficients disagree due to the very contribution from the 
evanescent operators. We have cross-checked our results using light quark masses as IR 
regulators, in which case there was no need to consider evanescent operators. 

The effect of the NLO QCD corrections on the branching ratios varies between and 
around 15% in a reasonable class of renormalization schemes. Higher order effects are ex- 
pected to be around 1%. The effect of the evanescent operator contribution in the matching 
turns out to be small. Thus, from the numerical point of view, our results are a confirmation 
of what has been found in ref. M and in the following phenomenological analyses [O . 
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